Maple 2025 Questions and Posts

These are Posts and Questions associated with the product, Maple 2025

I wanted to check that the input  has the pattern   symbol(symbol), which will match any of   y(x), or f(x) or A(B) and so on.

But using patmatch does not work. Using patmatch(h(z),y::symbol(x::symbol),'la');  or even patmatch(h(z),'y'::anything('x'::anything),'la'); all return false. I know I can do patmatch(h(z),func::function(name),'la'); and this returns true, but this matches h(z,r) and matches h(z,r,t) and matches h(z,r,t,u) and so on. 

I wanted to match only   SYMBOL(SYMBOL), i..e. one symbol followed by "(" followed by one symbol followed by closing ")"

For reference, this is what I am looking for 

I know I can use other ways in Maple to do this (may be typematch and and others). But wanted to see if patmatch works on this and why it is failing.

Can this be done using patmatch?

Maple 2025.1 unable to solve this ode. Sympy gives the following two solutions which Maples verifies are correct.

Any trick or option that can help dsolve find these solutions?
 

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

restart;

ode:=diff(y(x),x) = (1+cos(x)*sin(y(x)))*tan(y(x));

diff(y(x), x) = (1+cos(x)*sin(y(x)))*tan(y(x))

sol:=dsolve(ode);

sol_1:=y(x)=arcsin( 2*exp(x) / ( c__1 + sqrt(2)*exp(x) * sin(x+Pi/4) ) ) + Pi

y(x) = arcsin(2*exp(x)/(c__1+2^(1/2)*exp(x)*sin(x+(1/4)*Pi)))+Pi

odetest(sol_1,ode)

0

sol_2:=y(x)=arcsin( 2*exp(x) / ( c__1 - sqrt(2)*exp(x) * sin(x+Pi/4) ) ) ;

y(x) = arcsin(2*exp(x)/(c__1-2^(1/2)*exp(x)*sin(x+(1/4)*Pi)))

odetest(sol_2,ode)

0

 


 

Download How_to_find_solution_sept_20_2025.mw

update:

OK, found out how. Needed transformation u(x)=sin(y(x)). Maple probably did not have this in one of the things to try.

 

restart;

ode:=diff(y(x),x) = (1+cos(x)*sin(y(x)))*tan(y(x));
sol:=dsolve(ode);

diff(y(x), x) = (1+cos(x)*sin(y(x)))*tan(y(x))

tr:=y(x)=arcsin(u(x));
PDEtools:-dchange(tr,ode,[u(x)]):
dsolve(%);
sol:=y(x)=arcsin(rhs(%));
odetest(sol,ode)
 

y(x) = arcsin(u(x))

u(x) = -2/(-2*exp(-x)*c__1+sin(x)+cos(x))

y(x) = -arcsin(2/(-2*exp(-x)*c__1+sin(x)+cos(x)))

0


 

Download How_to_find_solution_sept_20_2025_V2.mw

 

 

Why when given IC for this ode, where the IC do not really makes much sense, so was not used. But the question is on the format of the output of the Maple dsolve. It gives solution as [{y(t) = c__1}]  instead of y(t) = c__1
 

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

restart;

ode:=diff(y(t),t)=0;
IC:=y(0)=t;
sol:=dsolve(ode)

diff(y(t), t) = 0

y(0) = t

y(t) = c__1

sol:=dsolve([ode,IC])

[{y(t) = c__1}]


 

Related question. Since Maple did not use the IC, should there have been warning message generated that IC was ignored?

 

Download strange_format_of_solution_sept_19_2025.mw

 

 

Maple documents appear to be an excellent way of documenting problems and procedures, and notetaking in general, but I have always found them laborious and frequently encountered problems with formatting. One may as well use Tex. This is the reason that I was struck by the ability of the GenerateDocument() procedure in the NaturalLanguage package to create Maple documents.

Is there any practical way to use NaturalLanguage for notetaking and documentation of procedures. By this I mean, the ability to input math and or explanatory text and request(direct) the AI to produce a maple document.

I have already plotted the Nusselt number as a line graph. Now, I would like to plot it as a bar graph instead of line graphs. Specifically, I want to create grouped bar charts that combine several parameters (such as Gr, Rb, N[t], N[r], M, Sc) in a single figure. Each group should correspond to one parameter, with bars representing distinct values (e.g., 0.1, 0.4, 0.6, 0.8). For reference, I have attached a sample figure from another work. Could anyone please help in this regard?

Group_bar_graphs_help.mw

Suppose I have a metric g, and I want to perform a conformal transformation g = exp(2Phi(X))*g, is there a straightforward way to do this for curvature quantitieies (Christoffel, Ricci Scalar etc)? I was able to do it rather easily for the Christoffel symbols, as seen below, but it required me making a substitution for each index pair. While this isn't horrible, it would be nice if there was a way to do it without this procedure.

** Edited to make it Phi(X) 

restart;

with(Physics):CompactDisplay(Phi(X)):

Phi(X)*`will now be displayed as`*Phi

(1)

Setup(metric=arbitrary):

_______________________________________________________

 

`Systems of spacetime coordinates are:`*{X = (x1, x2, x3, x4)}

 

`Default differentiation variables for d_, D_ and dAlembertian are:`*{X = (x1, x2, x3, x4)}

 

`The arbitrary metric in coordinates `*[x1, x2, x3, x4]

 

`Signature: `(`- - - +`)

 

_______________________________________________________

 

Physics:-g_[mu, nu] = Matrix(%id = 36893488152300428092)

(2)

Christoffel[~alpha,beta,nu,definition]

Physics:-Christoffel[`~alpha`, beta, nu] = (1/2)*Physics:-g_[`~alpha`, `~mu`]*(Physics:-d_[nu](Physics:-g_[mu, beta], [X])+Physics:-d_[beta](Physics:-g_[mu, nu], [X])-Physics:-d_[mu](Physics:-g_[beta, nu], [X]))

(3)

Simplify(eval((3),[g_[mu,beta]=exp(2*Phi(X))*g_[mu,beta],g_[mu,nu]=exp(2*Phi(X))*g_[mu,nu],g_[beta,nu]=exp(2*Phi(X))*g_[beta,nu]]))

Physics:-Christoffel[`~alpha`, beta, nu] = exp(2*Phi(X))*(Physics:-d_[nu](Phi(X), [X])*Physics:-g_[beta, `~alpha`]+Physics:-d_[beta](Phi(X), [X])*Physics:-g_[nu, `~alpha`]-Physics:-d_[`~alpha`](Phi(X), [X])*Physics:-g_[beta, nu]+Physics:-Christoffel[`~alpha`, beta, nu])

(4)

expand((4))

Physics:-Christoffel[`~alpha`, beta, nu] = (exp(Phi(X)))^2*Physics:-d_[nu](Phi(X), [X])*Physics:-g_[beta, `~alpha`]+(exp(Phi(X)))^2*Physics:-d_[beta](Phi(X), [X])*Physics:-g_[nu, `~alpha`]-(exp(Phi(X)))^2*Physics:-d_[`~alpha`](Phi(X), [X])*Physics:-g_[beta, nu]+(exp(Phi(X)))^2*Physics:-Christoffel[`~alpha`, beta, nu]

(5)
 

 

Any thoughts appreciated, thank you!

MyConformal.mw

Hello everyone!
I have had an issue for the past weeks, where it seems like Maple has a problem identifying the correct units and also sometimes having an issue with defining a variable.

In the picture below you can see I have defined rho, m and tried to Solve V. I get m^3, which is fine but I cannot change the unit in the right bar (see picture).

I even tried to just take square root of my V to see if I could then change units. It seems like it thinks I'm playing with weight.
I even had 2 teachers trying to help me find a solution, without luck.

I have tried executing the whole document and also only bits of the document without luck.

Does anyone know this problem and has a solution?

Thanks in advance! 

- Jacob

For the derivation of the results that I need for my research, I use advanced calculus, special functions etc.

Consider a derivation of the form LHS = RHS1 = ... RHSn.

Suppose that I have a LaTeX file describing this derivation, with fine details of the derivation in notes written using pencil on paper. 

To check this derivation for errors, such as a wrong sign, I could retype LHS and RHSn into Maple and check for equality using symbolic algebra.

To provide a more detailed check of this derivation, I could type LHS, RHS1, ..., RHSn into Maple and check the sequence equalities LHS = RHS1, ..., RHS(n-1) = RHSn using symbolic algebra.

However, this is a laborious and, itself, an error-prone process.

Is there some way of submitting my LaTeX file (or its compiled version) and/or an image of my handwritten notes for checking by Maple, line-by-line, using symbolic algebra?

These are current collections of Maple bugs before I lose track of them. I put them all in one post. Hopefully these can be fixed in Maple 2025.2. For each problem, I post separate worksheet, so there are few worksheets here.

This is all on Linux using 2025.1 and latest SupportTools and latest Physics.

1. Random crashes. This one is very strange. The crash happens randomly. You might need to try few times to see it or close the worksheet and reopen it.
 

restart;

Example . RANDOM CRASHES

 

restart;

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC],y(t)) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

[(5/2)*2^(1/2)*(-(1-sin(2*t)*sin(8)-cos(2*t)*cos(8))^(1/2)*(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2)-sin(2*t)*cos(8)+sin(8)*cos(2*t))/(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2), -10]

restart;

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC],y(t)) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

Error, (in anonymous procedure called from cos) too many levels of recursion

restart;

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC]) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

[-(5/2)*2^(1/2)*((1-sin(2*t)*sin(8)-cos(2*t)*cos(8))^(1/2)*(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2)-sin(8)*cos(2*t)+sin(2*t)*cos(8))/(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2), -10]

restart;

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC]) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

[-(5/2)*2^(1/2)*((1-sin(2*t)*sin(8)-cos(2*t)*cos(8))^(1/2)*(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2)-sin(8)*cos(2*t)+sin(2*t)*cos(8))/(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2), -10]

restart;

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC]) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

Error, (in signum) too many levels of recursion

restart;

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC],y(t)) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

Error, (in anonymous procedure called from cos) too many levels of recursion

restart

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC],y(t)) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

[(5/2)*2^(1/2)*(-(1-sin(2*t)*sin(8)-cos(2*t)*cos(8))^(1/2)*(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2)+sin(8)*cos(2*t)-sin(2*t)*cos(8))/(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2), -10]

restart;

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC],y(t)) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

[-(5/2)*2^(1/2)*((1-sin(2*t)*sin(8)-cos(2*t)*cos(8))^(1/2)*(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2)+sin(2*t)*cos(8)-sin(8)*cos(2*t))/(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2), -10]

restart;

sol:=y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2);
ode:=diff(y(t),t) = (25-y(t)^2)^(1/2);
IC:=y(4)=-5;
odetest(sol,[ode,IC]) assuming t>1;

y(t) = (-25*sin(t+arctan(sin(4)*cos(4)/(sin(4)^2-1)))^2+25)^(1/2)

diff(y(t), t) = (25-y(t)^2)^(1/2)

y(4) = -5

[-(5/2)*2^(1/2)*((1-sin(2*t)*sin(8)-cos(2*t)*cos(8))^(1/2)*(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2)+sin(2*t)*cos(8)-sin(8)*cos(2*t))/(1+sin(2*t)*sin(8)+cos(2*t)*cos(8))^(1/2), -10]

 

 

 

 

Download random_crashes_sept_8_2025.mw

 

2. collection of bugs from solve(identity) (another one related to solve(identity at end)

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1877 and is the same as the version installed in this computer, created 2025, July 11, 19:24 hours Pacific Time.`

 

Example 1

 

restart;

eq:=1/8*A^2*exp(2*theta*(B+I))+1/8*exp(2*theta*(B-I))*A^2-1/4*A^2*exp(2*B*theta)-1/4*exp(theta*(B-2*I))*A*B-1/4*exp(theta*(B+2*I))*A*B+1/2*A*B*exp(B*theta)+1/4*exp(theta*(B-2*I))*A*C+1/4*exp(theta*(B+2*I))*A*C-1/2*A*C*exp(B*theta)-1/4*I*exp(theta*(B-2*I))*A+1/4*I*exp(theta*(B+2*I))*A+1/4*C^2*cos(2*theta)-1/4*C^2-1/2*C*sin(2*theta)-1/2*cos(2*theta)-1=0:
the_vars:=[A, B, C]:
solve(identity(eq,theta),the_vars);

Error, (in gcd/doit) too many levels of recursion

 

Example 2

 

restart;

eq:=-x^(1/2)-1/2*x*A^2+A*B*sinh(B*x)-1/2*x*A^2*cosh(2*B*x)=0;
the_vars:=[A, B]:
solve(identity(eq,x),the_vars);

-x^(1/2)-(1/2)*x*A^2+A*B*sinh(B*x)-(1/2)*x*A^2*cosh(2*B*x) = 0

Error, (in gcd/doit) too many levels of recursion

 

 

Example 3

 

restart;

eq:=1 = X*(2*cos(X)*cos(x0)-X*sin(X)*cos(x0)-2*sin(X)*sin(x0)-X*cos(X)*sin(x0)-x0*sin(X)*cos(x0)-x0*cos(X)*sin(x0))*(2*Y*ln(Y+y0)+Y+2*y0*ln(Y+y0)+y0)/Y/(X*cos(X)*cos(x0)-X*sin(X)*sin(x0)+x0*cos(X)*cos(x0)-x0*sin(X)*sin(x0)+sin(X)*cos(x0)+cos(X)*sin(x0))/(2*ln(Y+y0)+2*Y/(Y+y0)+1+2*y0/(Y+y0));

1 = X*(2*cos(X)*cos(x0)-X*sin(X)*cos(x0)-2*sin(X)*sin(x0)-X*cos(X)*sin(x0)-x0*sin(X)*cos(x0)-x0*cos(X)*sin(x0))*(2*Y*ln(Y+y0)+Y+2*y0*ln(Y+y0)+y0)/(Y*(X*cos(X)*cos(x0)-X*sin(X)*sin(x0)+x0*cos(X)*cos(x0)-x0*sin(X)*sin(x0)+sin(X)*cos(x0)+cos(X)*sin(x0))*(2*ln(Y+y0)+2*Y/(Y+y0)+1+2*y0/(Y+y0)))

solve(identity(eq,X),[x0,y0]);

Error, (in signature) too many levels of recursion

solve(identity(eq,X),[x0,y0,Y]);

Error, (in signature) too many levels of recursion

 

 


 

Download collection_of_maple_internal_errors_sept_6_2025.mw

 

3. Adding Physics:-Setup(assumingusesAssume = true): make combine fail

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1877 and is the same as the version installed in this computer, created 2025, July 11, 19:24 hours Pacific Time.`

restart

Physics:-Setup(assumingusesAssume = true):

A:=1/6*ln(u^2+1)+1/3*arctan(u)+1/6*ln(u^2-3^(1/2)*u+1)-1/3*arctan(2*u-3^(1/2))+1/6*ln(u^2+3^(1/2)*u+1)-1/3*arctan(2*u+3^(1/2));
combine(A,ln) assuming real;

(1/6)*ln(u^2+1)+(1/3)*arctan(u)+(1/6)*ln(u^2-3^(1/2)*u+1)-(1/3)*arctan(2*u-3^(1/2))+(1/6)*ln(u^2+3^(1/2)*u+1)-(1/3)*arctan(2*u+3^(1/2))

Error, (in assuming) when calling 'is'. Received: 'invalid input: (u^2+1)^(1/6)*(u^2-3^(1/2)*u+1)^(1/6) <> 0'

Physics:-Setup(assumingusesAssume = false):

combine(A,ln) assuming real;

ln((u^2+1)^(1/6)*(u^2-3^(1/2)*u+1)^(1/6))+ln((u^2+3^(1/2)*u+1)^(1/6))+(1/3)*arctan(u)-(1/3)*arctan(2*u-3^(1/2))-(1/3)*arctan(2*u+3^(1/2))

 


 

Download adding_Phsyics_makes_combine_fail_sept_6_2025.mw

 

4. odetest internal error when adding assuming

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1877 and is the same as the version installed in this computer, created 2025, July 11, 19:24 hours Pacific Time.`

restart;

sol:=y(x) = 6*x/(3*x-2*LambertW(-3/2*exp(5/2*x+5/6*_C2)))+1/2*x+1/3;
ode:=x-2*y(x)-1+(3*x-6*y(x)+2)*diff(y(x),x) = 0;
odetest(sol,ode,y(x)) assuming positive;

y(x) = 6*x/(3*x-2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*_C2)))+(1/2)*x+1/3

x-2*y(x)-1+(3*x-6*y(x)+2)*(diff(y(x), x)) = 0

Error, (in depends) too many levels of recursion

odetest(sol,ode,y(x)); #removing positive it now works

-(40/3)*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))^4/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))+180*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))^3*x/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))-450*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))^2*x^2/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))+315*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))*x^3/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))-(40/3)*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))^3/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))-252*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))^2*x/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))+630*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))*x^2/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))+315*x^3/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))-432*x*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))/((-3*x+2*LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2)))^3*(1+LambertW(-(3/2)*exp((5/2)*x+(5/6)*c__2))))

 


 

Download internal_odetest_error_sept_6_2025.mw

 

5. solve(identity,..  gives internal error when one variable is missing

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1877 and is the same as the version installed in this computer, created 2025, July 11, 19:24 hours Pacific Time.`

restart;

eq:=-A^2*exp(2*B*x)+A*B*exp(B*x)-2*A*C*exp(B*x)-C^2-a*cos(b*x)^m*(A*exp(B*x)+C+1)=0;

-A^2*exp(2*B*x)+A*B*exp(B*x)-2*A*C*exp(B*x)-C^2-a*cos(b*x)^m*(A*exp(B*x)+C+1) = 0

the_vars:=[A, B, C,m]: #all variables are listed
solve(identity(eq,x),the_vars);

[[A = 0, B = B, C = -(1/2)*a-(1/2)*(a^2-4*a)^(1/2), m = 0], [A = 0, B = B, C = -(1/2)*a+(1/2)*(a^2-4*a)^(1/2), m = 0], [A = -C-(1/2)*a-(1/2)*(a^2-4*a)^(1/2), B = 0, C = C, m = 0], [A = -C-(1/2)*a+(1/2)*(a^2-4*a)^(1/2), B = 0, C = C, m = 0]]

the_vars:=[A, B, C]:   #forget to add m variable to list, now it gives internal error variables are listed
solve(identity(eq,x),the_vars);

Error, (in depends) too many levels of recursion

 


 

Download missing_variable_solve_sept_6_2025.mw

 

6. odesteps gives internal error (was question before, moved it to here, so all in one place)
 

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

restart;

ode:=x^2*diff(y(x),x$2)+(x^2-5*x)*diff(y(x),x)+(5-6*x)*y(x)=0; #22942.  

x^2*(diff(diff(y(x), x), x))+(x^2-5*x)*(diff(y(x), x))+(5-6*x)*y(x) = 0

sol:=dsolve(ode);

y(x) = c__1*x^5*(x+5)+c__2*x*(x^4*(x+5)*Ei(1, x)+(-x^4-4*x^3+3*x^2-4*x+6)*exp(-x))

Student:-ODEs:-ODESteps(ode)

Warning, cannot verify that the given particular solution, y(x) = 1+1/5*x, actually solves the corresponding homogeneous ODE, diff(diff(y(x),x),x)+1/x*(x-5)*diff(y(x),x)-(-5+6*x)/x^2*y(x) = 0

Error, (in Student:-ODEs:-ChangeVariables) the ODE, diff(diff(U(T),T),T) = 5*(T^2+6*T-5)/T^2/(5+T)*U(T)-diff(U(T),T)*(T^2+2*T-25)/T/(5+T), contains the undifferentiated dependent variable, U(T), but the transformation %3, does not

 


 

Download internal_error_ODESteps_sept_2_2025.mw

I'm not sure exactly when this started but maple on one of my linux machines won't print more than one row of a matrix. I get dots like this:

Note that interface(rtablesize) is [10,10]. This happens with both maple2024 and maple2025. 

On any of my other machines, maple works just fine. I have also tried deleting ~/.maple and ~/.maplesoft

Any idea why Maple dsolve can't find solution to this ode? From textbook

The strange thing, it solves if it asked for implicit solution. But the default, will give no solution.

Is this a defect? Should it not have returned the book solution automatically?   How is a user supposed to know the ode has a solution or not, if default call returns no solution?

restart;

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

Physics:-Version();

`The "Physics Updates" version in the MapleCloud is 1877 and is the same as the version installed in this computer, created 2025, July 11, 19:24 hours Pacific Time.`

restart;

ode:=v(x)*diff(v(x),x) = g;
ic:=v(x__0) = v__0;
sol:=dsolve([ode,ic]);

v(x)*(diff(v(x), x)) = g

v(x__0) = v__0

restart;

ode:=v(x)*diff(v(x),x) = g;
ic:=v(x__0) = v__0;
sol:=dsolve([ode,ic],'implicit');

v(x)*(diff(v(x), x)) = g

v(x__0) = v__0

-2*g*x+v(x)^2+2*g*x__0-v__0^2 = 0

#why did not default call return this?
PDEtools:-Solve(sol,v(x))

v(x) = (2*g*x-2*g*x__0+v__0^2)^(1/2), v(x) = -(2*g*x-2*g*x__0+v__0^2)^(1/2)

Download dsolve_gives_no_solution_sept_2_2025.mw

The_Bohrs_Model_-_MaplePrimes.mw

Look at the equation (11) in the Maple's document. I would like to force Maple to let the variable "r" inside the squared root so to get the equation (12). Any idea of doing that?  Thank you in advance for your help.

(2500iw/(1+5iw) )+(200iw/1-10iw)+5 rationalize and simplify

I was surprised that Maple can't solve this first order ode which is exact ode.

I solved by hand and Maple says my solution is correct.

Any one can find why Maple failed to solve this and if older versions can solve it? Also tried implicit option, but that did not help.

interface(version);

`Standard Worksheet Interface, Maple 2025.1, Linux, June 12 2025 Build ID 1932578`

SupportTools:-Version();

`The Customer Support Updates version in the MapleCloud is 29 and is the same as the version installed in this computer, created June 23, 2025, 10:25 hours Eastern Time.`

restart;

ode:=diff(y(x),x) = (2*sin(2*x)-tan(y(x)))/x/sec(y(x))^2;

diff(y(x), x) = (2*sin(2*x)-tan(y(x)))/(x*sec(y(x))^2)

sol:=dsolve(ode);

mysol:=cos(2*x)+x*tan(y(x))=c__1;

cos(2*x)+x*tan(y(x)) = c__1

odetest(mysol,ode);

0

 

 

Download maple_solving_exact_ode_august_25_2025.mw

As we can see, RealDomain:-solve gives an incorrect solution to the following system: 

restart;

sys := `~`[diff](sqrt(2*a^2-8*a+10)+sqrt(b^2-6*b+10)+sqrt(2*a^2-2*a*b+b^2), [a, b]):

RealDomain:-solve(`~`[`=`](sys, 0), {a, b})

{a = 5/3, b = 5/2}, {a = a, b = 2*a/(a-1)}

(1)

plot(eval(sys, {max(2*5^(1/2), (2*a^2-8*a+10)^(1/2)+2^(1/2)*((a^2-4*a+5)/(a-1)^2)^(1/2)+2^(1/2)*(a^2*(a^2-4*a+5)/(a-1)^2)^(1/2)), min(2*5^(1/2), (2*a^2-8*a+10)^(1/2)+2^(1/2)*((a^2-4*a+5)/(a-1)^2)^(1/2)+2^(1/2)*(a^2*(a^2-4*a+5)/(a-1)^2)^(1/2))}[-1]), a = -infinity .. infinity)

 

extrema(sqrt(2*a^2-8*a+10)+sqrt(b^2-6*b+10)+sqrt(2*a^2-2*a*b+b^2), {}, {a, b})

{max(2*5^(1/2), (2*a^2-8*a+10)^(1/2)+2^(1/2)*((a^2-4*a+5)/(a-1)^2)^(1/2)+2^(1/2)*(a^2*(a^2-4*a+5)/(a-1)^2)^(1/2)), min(2*5^(1/2), (2*a^2-8*a+10)^(1/2)+2^(1/2)*((a^2-4*a+5)/(a-1)^2)^(1/2)+2^(1/2)*(a^2*(a^2-4*a+5)/(a-1)^2)^(1/2))}

(2)

Download solve_returns_an_unsatisfiable_real_solution.mw

This appears to be a bug; is it possible to fix it? 
Text: 

sys := diff~(sqrt(2*a^2 - 8*a + 10) + sqrt(b^2 - 6*b + 10) + sqrt(2*a^2 - 2*a*b + b^2), [a, b]):
RealDomain:-solve(sys =~ 0, {a, b});
1 2 3 4 5 6 7 Last Page 1 of 10